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Abstract 

In this work, the quantum confinement effect is proposed as the cause of the dis¬ 
placement of the vibrational spectrum of molecular groups that involve hydrogen bonds. 
In this approach the hydrogen bond imposes a space barrier to hydrogen and constrains 
its oscillatory motion. We studied the vibrational transitions through the Morse po¬ 
tential, for the NH and OH molecular groups inside macromolecules in situation of con¬ 
finement (when hydrogen bonding is formed) and non-confinement (when there is no 
hydrogen bonding). The energies were obtained through the variational method with 
the trial wave functions obtained from Supersymmetric Quantum Mechanics (SQM) 
formalism. The results indicate that it is possible to distinguish the emission peaks 
related to the existence of the hydrogen bonds. These analytical results were satisfac¬ 
torily compared with experimental results obtained from infrared spectroscopy. 


* Institute de Biociencias, Letras e Ciencias Exatas, IBILCE-UNESP, Sao Jose do Rio Preto, SP, Brazil 
Hnstituto de Biociencias, Letras e Ciencias Exatas, IBILCE-UNESP, Sao Jose do Rio Preto, SP, Brazil 
Raculdade de Tecnologia de Sao Paulo, FATEC/SP-CEETPS-UNESP, Sao Paulo, SP, Brazil 


1 



INTRODUCTION 


Hydrogen bond is an important interaction present in different systems in nature; the most 
famous examples of its occurrence are in water molecules and the DNA double helix, ^ and 
references therein. Particularly, biological macromolecules as proteins, nucleic acids and car¬ 
boxylic acids have the molecular groups NH and OH in their structures. These groups induce 
the hydrogen bonds which are essential to stabilize the structures and play an important role 
in the behavior of these macromolecules. Several other systems like clathrates and liquid 
crystals present hydrogen bonds in their structures^. 

Few years ago a modern dehnition of the hydrogen bond following a list of criteria was 
recommended, after a survey on recent results of the subject,^. According to this dehnition 
the hydrogen bond is an attractive interaction between a hydrogen atom from a molecule or a 
molecular fragment X-H in which X is more electronegative than H, and an atom or a group 
of atoms in the same or a different molecule, in which there is evidence of bond formation. 
The hydrogen bond can be classihed as a donor-acceptor interaction. The notation of the 
bond is X-H- ■ ■ Y, where the dots represent the bond, X-H is the hydrogen bond donor and 
Y is the acceptor anion. The covalent bond is represented by X-H and the hydrogen bond is 
represented by H- ■ ■ Y. The high electro-negativity of atom X compared with the hydrogen 
atom in the X-H chemical bond makes the hydrogen to share part of its electronic density. 

The different spectroscopy techniques, such as, Raman, infrared (IR), Terahertz (THz) 
and Nuclear Magnetic Resonance (NMR) spectroscopy, provide the identihcation of the 
vibrational transitions of the molecular groups involved, since they are sensitive to changes 
in the spectrum. Therefore, it is possible to know if a specihc molecular group participate 
of the hydrogen bond, because the vibrational spectrum changes. 

In this work it is proposed that quantum conhnement effect is the main cause of de¬ 
localization of the vibrational spectrum of NH and OH groups when hydrogen bonds are 
formed. 

Quantum conhnement can be dehned as the restricted space where the particles can 
move and it is generally introduced by a potential barrier created, for instance, due to elec¬ 
trostatic potentials, the presence of an interface between different materials or the presence 
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of a surface. The simplest case of confinement consists of a particle subject to an infinite 
square well potential treated in introductory quantum mechanics,®. In this confined regime, 
the motion of a particle is restricted to a well defined region of space and there is the ap¬ 
pearance of discrete energy levels due to quantization effects. They appear once the size 
of the barrier is of the same magnitude as the de Broglie wavelength of the particle wave 
function. Conhned systems have interesting applications; in material science the quantum 
conhnement effect is well studied theoretically and experimentally in semiconductor physics, 
mainly in the study of quantum dots,®,®. If the size of the material is reduced, the elec¬ 
tronic and optical properties are drastically changed and can be controlled; this is one of the 
main subjects of nanostructure science,®. There are other cases that may also be studied 
by confinement as, for instance, a system subjected to high pressure,®. Recent studies focus 
on theoretical as well as experimental investigations of hydrogen/halogen bonds in different 
conhning environments,®-^. 

In the approach used in this work the confinement is introduced by a potential barrier. 
In the formation of the hydrogen bond, the acceptor atom of the hydrogen bond introduces 
a barrier which prevents the free oscillation of the hydrogen atom attached to the donor 
atom. In this situation, the hydrogen is confined between two atoms and the energy levels 
are distorted by this conhnement, i. e., there are vibrational frequency shifts. Vibrational 
frequency shifts can also be induced by high pressure,^ -1^. 

Thus, the main interest of this paper is to study the hydrogen bonds present in biological 
macromolecules, particularly, to determine the vibrational transitions related to OH and 
NH groups. To overcome this proposal we used an analytical methodology based on the 
variational method associated to Supersymmetric Quantum Mechanics, (SQM). 

The conhned quantum mechanical problem is settled by considering the Schrddinger 
equation depending on a potential with a potential barrier. Considering that the hydrogen 
bonds are highly directional, the vibrations of the groups studied are better described by 
the one-dimensional Morse potential. If not by the potential barrier (without conhnement) 
this potential presents exact analytical solution,!^. However, by considering the conhned 
case (with conhnement) the problem is not exactly solvable and therefore an approximation 
method has to be introduced. At this point we use the formalism of SQM to calculate the 
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spectrum of wave functions and energies and use the variational method to treat the conhned 
case, with a potential barrier. The transition energy between two adjacent states is calculated 
and the values obtained are compared with experimental spectroscopy data to validate the 
proposed hypothesis. The formalism of SQM coupled to the variational method has been a 
successful methodology to study many problems, in particular, the conhned systems 

Other theoretical models based on the use of the Morse potential have been proposed to 
describe high frequency and low frequency vibration modes of the hydrogen bonds,®-l25l 
the approach proposed here, unlike these models, there is no need either to use any coupling 
theory nor to enter new parameters in the Morse potential. The data used here refer only 
to the covalent bonding and to the length of the hydrogen bond, which reduces signihcantly 
the number of parameters. 

MODEL 

A schematic model of a hydrogen bond between atoms X and Y is shown in Figure 1. 

In a hydrogen bonding, the hydrogen nuclei may oscillate within a region determined between 
a minimum distance (xmin) and a maximum distance (xmax)- The value of Xmin is given by 
the covalent radius (xc) of atom X and Xmax is given by the length of the hydrogen bond 
(lx--y) subtracted by the van der Waals radius of the acceptor atom, i.e., 

Xmax ^X-'-Y X^ix (1) 

The origin of the coordinate system is located at the center of the atom X. Any distance 
used is computed from this point and it is considered center-to-center. The covalent radius 
is adopted by the condition of impenetrability of the hydrogen atom donor; the employment 
of the van der Waals radius is justihed because it is a measure of the volume excluded by the 
acceptor atom, i.e., the hydrogen cannot penetrate this region. The atoms are considered 
as hard spheres and the hydrogen is treated as a material point, since the electron cloud is 
shared with the donor and the acceptor atoms. 

As mentioned before, the vibrational energy of the system is described through the con¬ 
hned one-dimensional Morse potential. This potential has been widely used to describe 
oscillations of diatomic molecules,®-^. The energy absorbed/emitted is obtained from the 
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calculation of the energies of the ground state (n = 0) (for NH and OH), the hrst excited 
state (n = 1) (for NH and OH) and the second excited state (n = 2) (only for OH). The 
energy of the absorbed/emitted photon is obtained from the difference between two adja¬ 
cent levels. The energies of the conhned case are evaluated by using the variational method 
associated to the formalism of SQM. 


SUPERSYMMETRIC QUANTUM MECHANICS, SQM 


SQM is an algebraic method commonly used to exploit different aspects of non-relativistic 
quantum mechanical systems. It is particularly efficient to solve exactly all the quantum 
potential problems by means of the Hamiltonian factorization, thus providing the entire 
spectrum of wave functions with their respective energies together with a hierarchy of Hamil¬ 


tonians, all related by the supersymmetric algebra,^. On the other hand, when the potential 


is non-exactly solvable, such as in a conhned system, an approximation technique is needed 


and the variational method has already appeared as fully appropriate,^. 


Thus, using the superalgebra a given Hamiltonian can be factorized in terms of bosonic 
operators. In h = c = 1 units, the Schrddinger equation can be rewritten as 



( 2 ) 


where is the lowest eigenvalue with ground state wave function satisfying 


H-t[,^^(x) = 0 


(3) 


and with the bosonic operators dehned as 



(4) 


where the superpotential Wi{x) satishes the Riccati equation 



(5) 


and the ground state wave function has the following form 



( 6 ) 
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Thus a whole hierarchy of Hamiltonians can be constructed, with simple relations con¬ 
necting the eigenvalues and eigenfunctions of the n-members, 


//„= a:A- + 4">, 

(7) 


(8) 

= yi+A+ Ei'i = 

(9) 


In equation (|^, the upper index between parentheses refers to the hierarchy Hamiltonian 
member and the lower index to the level within the hierarchy. 


MATHEMATICAL FORMALISM 


The one-dimensional Morse potential^, is given by 




( 10 ) 


where De is the dissociation energy of molecule, (3 is the parameter related to the width of 
potential well and Xeq is the inter-nuclear equilibrium distance. The (3 parameter may be 
determined by treating the molecule as a frequency oscillator v and through experimental 
data. By expanding the Morse potential in a Taylor series around the equilibrium position 
to the second order we hnd that 


/^ = 



1 

2Z/ . 


( 11 ) 


The Schrddinger equation for a molecule subject to the Morse potential, according to the 
model, is written in one-dimensional Cartesian coordinates with the change of variables, 
y = f3x. After multiplying both sides by we obtain the Schrddinger equation 


dy^ 


+ A'^{e y-i')-2e ^"''^)T(?/) 


e^{y) 


( 12 ) 


where 


and 




2yDe 

(3^ A? 


2yE 


(13) 


(14) 
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where /i is the reduced mass of the system. The Hamiltonian operator in equation (12) is 


given by 

Hm = -^+ A2(e-2(3/-3/ed _ (15) 

dy^ 

Since the Schrddinger equation is exactly solvable, the process of factorization can be com¬ 
pletely made and the whole hierarchy can be evaluated. The result is,^ 

14 + 1 ( 2 /) = K^{e-^^y-y^o^ - 2e-^y-y^'>'^) + 


277 I 1 

n;„+i(i/) = + (A-^) 

4""“'=-(A - 44)^- 


( 16 ) 


In addition we can calculate exactly the wave functions. The ground state wave function is 
evaluated by using equation (§ and is given by 

Tq^^(i/) oc exp (^Ke~^y~y’"^'^ — y{K — -)^ . (17) 

This wave function is then used to evaluate, through the SQM algebra, the wave function 
for the excited states by using equation (9); the hrst excited state is given by 

oc —2 — A -|- l) exp — y{A — -)^ (18) 

and the second excited state is given by 

T^^^(i/) oc - SA^e-'^^y-y^’’^ + 12Ae-(2^-^^«) - A + dA^ - IdA + 12) . (19) 

exp {^—Ae~^y~y^'^^ — y{A — |)) . 


The wave functions given in (17), (18) and (19) are obtained to the exact system, non- 
conhned or with no hydrogen bond. We refer to them as '^n{non-confined)- For fhe confined 
system a potential barrier is introduced meaning that the wave functions must be zero 
outside the barrier. As the system is no longer exact, it is necessary to use an approximative 
method, which we choose to be the variational method. Thus, based on previous results,® 
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and^^, the eigenfunctions are multiplied by the term {ymax — y){y — ymin), i-e., the wave 
functions become zero at the edges of the conhnement region. For the ground state, this 
approach is equivalent to introduce two additional terms in the superpotential Wi, which 
corresponds to an inhnite barrier in the conhnement border, Vi(ymax) = Vi{ymin) = oo. 

Adopting this procedure, the wave functions are dehned for the region of space in which 
we are interested and, consequently, the probability of hnding the hydrogen outside the 
conhnement region must be zero. Thus, the trial wave functions to be used in the variational 
method are given by 


(confined) (non—confined) -(.ymax y')(,y ymin) 


( 20 ) 


and we will only need the cases of n = 0,1,2. 

The energies for the conhned case are obtained by solving 


En = 


I 


KHM'^ndy 

ymin 

^ ymin 


( 21 ) 


where Hm is, the Hamiltonian of the system given by equation (15) and are the trial wave 


functions given by equation (20), where n(non-confined) ^^^e given by (17), (18) and (19). The 


integration limits are the covalent radius of the atom for which the hydrogen is chemically 
bound {xmin) and the van der Waals radius complementary of hydrogen bond (xmax)- The 
results of energy are presented in cm~^ units (inverse wavelength) to facilitate comparison 
with the spectroscopic results. 


RESULTS 

The hydrogen bonds studied involve the NH and OH groups. Tables 1 and 2 below list the 
necessary parameters, obtained experimentally,^^. 


Table 1. Parameters used in the NH and OH groups energy evaluations,^. 








Equilibrium distance (xeq) 

Harmonic vibrational wavenumber (ca) 
Energy dissociation (De) 

Reduced mass (/r) 


NH 

OH 

1.0362 A 

0.9697 A 

3282cm“^ 

3738cm“^ 

3.47 eV 

4.39 eV 

1.5614.10-27fcc/ 

1.5746.10-27 


Table 2. Covalent and van der Waals radius used for the nitrogen and oxygen atoms, 


N O 


Covalent radius (xc) 0.71 A 0.66 A 

van der Waals radius (xv^) 1.55 A 1.52 A 


Results for NH group 

Non-confined case 


The f3 parameter in Morse potential is determined using the values given in Table 1 through 
the relation (|TT|) and for NH group it is given by 2.3168.10^°m“^. The value of the constant 


A for NH group, given by equation (13), is 17.0534. 

In the situation free of any conhnement, the calculation of the transition from the hrst 


excited state to the ground state is made according to equation (16). In this case, in spec¬ 
troscopic units, this value corresponds to AEi^q = Ei — Eq = 3089cm“^ and it should be 
checked by means of spectroscopy. By IR spectroscopy with polypeptides, an absorption 
band in 3090cm“^ is observed and it has been described as a vibration frequency of the NH 
molecular group,®'. This value is very close to the calculated and can be identihed, according 
to the proposed model, as vibration of NH without formation of hydrogen bonding. 
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Confined case 


When the NH group form molecular hydrogen bonding, the vibration frequency is shifted. 
An interesting system where this effect can be observed is in proteins. These may be in 
the form of a-helix, /3-sheet or random. In this helicoidal format the hydrogen bonding is 
between the amine group of an amino acid and the oxygen of the carboxyl group of another 
amino acid. Thus, proteins consist in a good system to verify the quantum conhnement effect 
in hydrogen bonding. It is found in the literature that proteins in a-helix format present hy¬ 
drogen bond lengths around 2.80A,I^. Considering the uncertainty of experimental results, 
the values of 2.75A, 2.80A and 2.85A were used as the maximum distance of the hydrogen 
oscillation, according to the equation Q. Solving equation (21) numerically, the results for 
the emission/absorption for proteins in a-helix are indicated below in Table 3. 


Table 3. Energy levels of proteins with hydrogen bonding, for different values of Xmax in 
(cm“^) units. 


^max 




2.75 A 

-26105 

-22643 

3461 

2.80A 

-26245 

-22835 

3410 

2.85A 

-26306 

-23010 

3296 


The hemoglobin, a metalloprotein present in the red blood cells of all vertebrates, is 
constituted of four a-helix chains. IR spectroscopy realized to verify the alterations caused by 
/3-thalassemia disorder, a reduction or absence synthesis of the beta chains of hemoglobin,!^, 
showed peaks around 3050cm“^ and 3290cm“^, that can be identihed as vibrations of the NH 
group in the situations free of conhnement and conhned, or with or without hydrogen bond, 
respectively. By considering 2.85A as the maximum distance of the hydrogen oscillation, 
the comparison of our theoretical result for the conhned case (3296cm“^) with these results 
shows an error of 0.18%. Thus we expect that other proteins in a-helix format present 
absorption peaks around 3296cm“^. In fact, other proteins in the a-helix format present 
transitions verihed by IR spectroscopy that corroborate with our results, as shown below. 

Keratin, for instance, is a protein found in abundance in life, in the human epidermis, hair. 
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nails, feathers, porcupine quills, etc. IR spectroscopy of virgin human hair (i.e., no chemical 
intervention) and horse hair, using lithium bromide and hydrochloric acid as reagent,^, 
showed the keratin absorption bands of hair with and without cuticle. Similar studies were 
carried out with a dried sample of goose feather followed by a soaked sample using heavy 
water,!^. The IR spectroscopy studies with some peptides, polyamides and hbrous proteins 
used swan feather wing and horse hair as a source of keratin,!^. 

Other protein having a-helix structure is rhodopsin, which can be found in the eye retina. 
IR spectroscopy with bovine rhodopsin Langmuir-Blodgett £lms,l^, showed an absorption 
peak in 3295cm“^; the theoretical result presents 0.03% error relative to this experimental 
result. 

The sodium-potassium pump (Na-f/K-l—ATPase) present in cells is another protein which 
contains several a-helices,^. The IR spectroscopy carried out in an aqueous medium of three 
different ionic compositions presented absorption peaks; in all types of solvent there is a peak 
located in 3285cm“^; the theoretical result presents 0.33% error relative to this experimental 
result. 

Table 4 summarizes the results discussed above. 

Table 4. Results for absorption/emission values = Ei — Eq, in {cm~^) units for the 

non-conhned and conhned cases with deviations relative to the experimental result. 
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Non-Confined Confined Deviation from 

(no hydrogen (with hydrogen experimental 



bond) 

bond) 

(%) 

Theoretical 

3089 

3296 


Hemoglobin with /3-thalassemia disorder,!^ 

3050 

3290 

0.18 

Normal Hair,l23l 

3066 

3300 

0.12 

Hair without cuticle,^ 

3062 

3295 

0.03 

Horse Hair,l23! 

- 

3293 

0.09 

Goose feather,l231 

3100 

3300 

0.12 

Swan Feather Wing,l^ 

3070 

3290 

0.18 

Horse Hair,l351 

3080 

3295 

0.03 

Rhodopsin (Retina eye),®' 

- 

3295 

0.03 

Sodium-Potassium Pump,®' 

- 

3285 

0.33 


Results for OH group 


Non-confined case 


In the case of OH gronp, the transition from the second to the hrst excited state (Aii^ 2 -?-i = 
E 2 — El) is also stndied. According to the formalism described, the (3 parameter, eqna- 


tion (11), for the OH gronp is 2.3558.10^^m ^ and the A constant valne is 18.9424. In the 
sitnation free of any type of conhnement, i.e. with no hydrogen bond, the transitions be¬ 


tween adjacent states were evalnated by using equation (16); the results are shown in Table 5. 


Table 5. Eigenvalues of energy to OH obtained by the model described in cm ^ units. 


^0 

El 

E 2 

^Ei^q 

Ai72_^i 

-33560 

-30019 

-26676 

3541 

3343 
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Confined case 


Carboxylic acid dimers are held by hydrogen bonds between the hydroxyl. By means of 
nentron diffraction, differently to the X-ray diffraction, one can obtain the distances between 
the hydrogen atom and the other atoms that compose the network,^. Figure 2 shows the 
distance between the hydrogen, linked to the donor atom, and the acceptor atom of the 
bond. Thus, to hnd the maximum distance of oscillation is necessary to sum the equilibrium 
distance of the OH bond to the length of hydrogen bond, i. e., the maximum distance is 
given by 

^max ^eq T ^H—O ^vw • 

Three types of carboxylic acids in the form of dimers, are studied in this work, the lengths 
of the hydrogen bonds are reported in Table 6. 


Table 6. Length of hydrogen bonds in dimers of different acids,^. 



Trihuoroacetic acid 

Acetic acid 

Formic acid 

Bond length ( Ih-o) 

1.658 A 

1.642 A 

1.660A 


Using the values of Table 6 to hnd the maximum distance of oscillation of hydrogen 


according to the equation (22) and solving equation (21) with wave functions given by (20) 


and (17)-(19) we hnd the value of emission/absorption energies for carboxylic acid dimers, 
shown in Table 7. 


Table 7. Energy levels of diherent systems, in [cm units. 



Eo 

E, 

E2 

Ai7i_j.o 

AF'2_s.i 

Trihuoroacetic acid 

-32953 

-29282 

-26270 

3671 

2944 

Acetic acid 

-32767 

-29342 

-26270 

3425 

3072 

Formic acid 

-32973 

-29277 

-26342 

3696 

2935 


The trihuoroacetic acid was analyzed in neon and argon matrices, in its monomeric 
form, in cyclic dimers and also in the deuterated form iCF^COOD) by means of IR spec- 
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troscopy,!^. In neon matrix the dimer of trifluoroacetic acid presented emission/absorption 
of 2922cm~^, described as the vibration of the OH molecular group, because when the acid 
is deuterated this peak disappears. The theoretical value presents 0.75% of error relative to 
this experimental value, considering the transition between the second excited state and the 
hrst excited state, AE 2 ^i = E 2 — Ei. 

Acetic acid was also analyzed by IR spectroscopy,^. There were treated several dimers 
of carboxylic acids. The values of 3067cm“^ in gas phase, 3100cm“^ carbon tetrachloride 
solution {CCI4) and 3080cm“^ in liquid phase agree with a transition AE 2 ^i (Table 7) of 
the conhned case, with theoretical errors smaller than 0.9%, obtained for the case where the 
acid was in CCl^ solution. 

IR spectroscopy was also used to analyse the formic acid photocatalytic decomposition 
process,®'. When it is analyzed on the catalyst surface of Ti02, without any additive, an 
absorption peak near 3698cm“^ can be noted, compatible to the transition AEi^q (Table 
7). This peak is attributed by the authors to the OH group, which during the reaction is 
displaced to values below 3668cm“^. This can be explained by the dissociation of the OH 
group during decomposition, preventing the formation of hydrogen bond and weakening the 
bonds still existing. 

The IF spectroscopy of photocatalytic decomposition of formic acid was analyzed in®l, 
showing transitions identihed as the vibration of the OH group, which agrees with the tran¬ 
sition of AEi^q (Table 7). Table 8 summarizes the results discussed above. 

Table 8. Results for absorption/emission values AEi^q = Ei — Eq, in {cm~^) units for the 
non-conhned and conhned cases with deviations relative to the experimental result. 
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AF'i-s.o 

Ai72-s>l 

Deviation from 

experimental (%) 

Triffuoroacetic acid,^ 

- 

2922 

0.75 

Theoretical 

3671 

2944 


Acetic acid,^ 




Gas phase 

- 

3067 

0.16 

In CCli solution 

- 

3100 

0.90 

Liquid phase 

- 

3080 

0.26 

Theoretical 

3425 

3072 


Formic acid,®' 

3698 

- 

0.05 

Formic acid,®' 

3697 

- 

0.03 

Theoretical 

3696 

2935 



CONCLUSIONS 

In this work the quantum conhnement effect in hydrogen bonds explains the shift in the 
emission spectrum of NH and OH groups. In all the cases mentioned the distortion in the 
spectrum when the bond is formed is considerable and measurable, and is easily distinguished 
from the non-confined case, with no hydrogen bond. This shift occurs because the hydrogen 
bond imposes a barrier to the hydrogen and prevents it from having freedom in its oscillatory 
motion. 

Quantitative results obtained for the vibrational transitions of the groups mentioned 
agree with the experimental results reported in the literature, reenforcing the applicability 
of the proposal and the reliability of this quantum mechanical methodology. The errors 
between the theoretical values entered by the variational method and the experimental data 
analyzed are small, the largest percentage error obtained in the calculations is smaller than 
1 . 0 %. 

The test wave functions used in the confined case, obtained with the aid of the SQM 
formalism, are appropriate to describe the phenomenon. As a final result of the construction 
of these wave functions, the polynomial terms are multiplied by the wave functions of the free 
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system, with no confinement. This term of confinement is sufficient to obtain the numerical 
results presented. 

As a final remark, the treatment proposed here can be extended to the study of other 
molecular groups, as FH and CH, involved in hydrogen bonds. 
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length of the hydrogen bond (/) 




Figure 1: Schematic one-dimensional model, where Xc is the covalent radius of the donor 
atom X, is the van der Waals radius of the acceptor atom Y and x^q is the equilibrium 
distance of covalent bond. 
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Figure 2: Schematic one-dimensional model for hydrogen bond for the group OH, remarking 
the length of hydrogen bond,*^. 
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